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Abstract. This paper is devoted to a non-perturbative renormalization group 
(NPRG) analysis of Model A, which stands as a paradigm for the study of 
critical dynamics. The NPRG formalism has appeared as a valuable theoretical 
tool to investigate non-equilibrium critical phenomena, yet the simplest — and 
nontrivial — models for critical dynamics have never been studied using NPRG 
techniques. In this paper we focus on Model A taking this opportunity to 
provide a pedagological introduction to NPRG methods for dynamical problems 
in statistical physics. The dynamical exponent z is computed in d = 3 and d = 2 
and is found in close agreement with results from other methods. 



1. Introduction 

The understanding of non-equilibrium critical phenomena stands as one of the major 
challenges of statistical physics. Systems far from thermal equilibrium arc omnipresent 
in nature: slow relaxation or external driving forces tend to prevent real systems from 
ever reaching their equilibrium distribution. Behavior out-of-equilibrium is far richer 
than at equilibrium, and many intriguing scaling phenomena, such as self-organized 
criticality (emergence of scaling without fine-tuning of a control parameter) pQ, or 
phase transitions between non-equilibrium stationary states 0] , have been observed 
for long. However, despite the considerable achievements of equilibrium statistical 
physics, the theoretical comprehension of non-equilibrium critical phenomena remains 
much poorer. The renormalization group (RG), which has appeared as a cornerstone to 
explain universality in equilibrium continuous phase transitions, has also allowed some 
breakthroughs out-of-equilibrium 0] . Nonetheless, many non-equilibrium phenomena 
remain out of range of perturbative approaches because of large coupling constants 
or because the interesting dimensions lie far from the critical one. Further theoretical 
progress out-of-equilibrium is hindered by the lack of analytical tools to handle the 
corresponding models. 

Recently, a novel approach - namely the non-perturbative renormalization group 
(NPRG) - has been proposed to investigate (non-equilibrium) reaction-diffusion 
processes [SJ. It has allowed to overcome the perturbative limitations and to gain 
physical insights into models such as branching and annihilating random walks - 
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which are reviewed in [S]. For instance, these studies have captured non-perturbative 
effects that essentially determine the phase diagram of some systems 00 
and have unveiled a genuinely non-perturbative fixed point governing the phase 
transition belonging to the so-called parity-conserving universality class [TD1 . A 
valuable advantage of this method is that it gives a unified description of a model: the 
very same equations enable one to probe any dimensions or coupling regimes, including 
non-universal features. Hence, the NPRG appears as a powerful tool to investigate 
non-equilibrium systems. 

However most readers are still largely unfamiliar with these techniques, though 
they have been introduced more than a decade ago for systems in equilibrium |11M12| . 
The aim of this contribution is to give a brief but pedagological introduction to the 
NPRG methods for non-equilibrium systems. For this purpose, Model A stands as 
one of the simplest - yet far from being trivial - dynamical models and it has never 
been studied within the NPRG framework so far, which we remedy with the present 
work. We adopt a practical viewpoint and put a particular emphasis on the discussion 
of the device of a (non-perturbative) approximation scheme. 

The paper is organized as follows. In Section [5] we briefly recall the definition of 
Model A and its critical properties. In Sectional the principles and the construction of 
the NPRG are outlined for generic non-equilibrium systems and specialized in Section^ 
to Model A. The corresponding non-perturbative flow equations are derived in Section 
[5] and their numerical integration is dealt with in Section [5] The results are eventually 
discussed in Section and followed by a brief summary in Section |H1 



2. Model A and critical dynamics 

The purely dissipative relaxation of a non-conserved field (f> can be described by the 
Langevin equation: 

dt cP(x,t)=-D-^^- + r ] (x,t) (1) 

where D denotes a constant and uniform relaxation rate and TL the usual Landau- 
Ginzburg- Wilson Hamiltonian. On approaching a critical point, the relaxation time 
of the order parameter starts diverging, which reflects the critical slowing down of the 
dynamics. The Langevin equation is a 'mesoscopic' description of the system which 
exploits the associated decoupling of time scales: the order parameter, represented by 
the (coarse-grained) field <j>, relaxes much slower than all the other microscopic degrees 
of freedom, which can hence be modeled by a stochastic Gaussian noise 77 with zero 
mean and variance 

(n(x,t)n(x',t')) = 2D k B T S d {x - x') S(t - t'). (2) 

The strength of the noise is fixed by the Einstein relation which ensures that the 
system acquires its equilibrium distribution at long time. We here focus on the case 
of a scalar order parameter with Ising symmetry, described by the Hamiltonian 

nm = J d d x Qiv^)] 2 + u(<f>) - h(x)<Kx)J , u(4>) = r -<j> 2 + ^cf> 4 . (3) 

The Langevin equation Q corresponds to a Glauber dynamics for the Ising spin 
and defines Model A in the classification by Halperin and Hohenberg |13|. Besides 
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the equilibrium critical exponents v and r?, the critical dynamics is characterized 
by the dynamical exponent z that describes the divergence of the relaxation time 
T r~j £ z ~ |T — T c \~ zl/ near the temperature T c of the second-order phase transition. 



Time Reversal Symmetry 

In the long-time limit after the initial perturbation, the system is expected to become 
time translational invariant (TTI) and the time-reversal symmetry (TRS) to hold. 
TRS then yields the fluctuation-dissipation theorem (FDT) which linearly relates the 
two-point correlation function C(x — x',t — t') — (<j)(x, t)<p(x', t')) with the response 
function R(x — x',t — t') = S((f>(x,t)) / '8h(x' ,t') following: 

R(x -x\t- t') = -6{t - t')-d t C(x -x',t- t'). (4) 

In the early stages of the relaxation process, the system is generally not TTI such 
that R and C may depend on both t and t' and FDT may not hold. The fluctuation- 
dissipation ratio TR/dtC becomes of particular interest to characterize the violation of 
FDT and the associated ageing phenomena (see m] for a recent review and references 
therein). For our non-perturbative study, we rather focus on the stationary dynamics 
where the system satisfies TRS. 



Field Theory 

Any RG treatment starts out from a field theory. Upon introducing a Martin-Siggia- 
Rose response field (j)(x,t) |15| . one can average over the Langevin noise 77 and cast 
the stochastic equation Q into a dynamic functional |16l I17| : 



S\6,$\ = / d d xdt 



8 t 6 + D § m 



I (5) 

(where fcgT has been set to unity). Correlation and response functions can then be 
expressed as functional averages with the weight exp(— S[cj>, </>]). 

In this equivalent field theoretical formulation, TRS can be conveniently expressed 
as an invariance of the action JSJ under a specific field transformation, as stressed in 
18 . This transformation writes 

Indeed, one can straigthforwardly check that after performing a time inversion t — > —t 
in J5j which switches the sign of the kinetic term 4>dt4>, the field transformation © 
yields additional contributions from the latter term and from the noise term D(j) 2 that 
cancel out. Besides, the transformation of the Hamiltonian part under @ produces an 
additional term oc dt<p STt/ '8<fi, which vanishes upon time integration in the stationary 
regime. We shall rely in the following on this simple expression of TRS to ensure that 
this invariance is preserved within the non-perturbative formulation. 



3. The NPRG formalism in non-equilibrium statistical physics 

The NPRG formalism relies on Wilson's RG idea [HJ\, which consists in building a 
sequence of scale-dependent effective Hamiltonians, that interpolate smoothly between 
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the short-distance physics at the (microscopic) scale k = A and the long-distance 
physics at the scale k — 0, through progressively averaging over fluctuations. Rather 
than expressing — as in the original Wilsonian formulation — the flow of effective 
Hamiltonians for the slow modes, one can work out the flow of effective 'Gibbs free 
energies' T k for the rapid ones, following [301 El- T/c thus only includes fluctuation 
modes with momenta \q\ > k. At the scale k = A, no fluctuation is yet taken into 
account and Ta coincides with the microscopic action S [23j, while at k = 0, all 
fluctuations are integrated out and To is the analogue of the Gibbs free energy T at 
thermal equilibrium, in that it encompasses the long-distance and long-time properties 
of the system. Hence, to construct T^, one needs at a given scale k to suppress the 
slow modes with momentum \q\ < k. This is achieved by adding to the original action 
JSJl a scale-dependent term [22 El 122] which is quadratic in the fields (so as to affect 
the propagator of the corresponding modes): 



AS k [<j>,4>] = l [ ^(x.tj.fet)]^^ 2 ,^)*^^),^,*)], (7) 

Z -lx,t 



where R k is a symmetric 2x2 matrix of elements (i,j — 1,2). These elements (so- 
called cutoff functions) will be specified in the following, but their general properties 
can be already stressed. In order to achieve the renormalization procedure outlined 
above, these cutoff functions must behave at fixed k as R k J ~ k 2 (in Fourier space) 
for small momenta |g| < k — so that the slow fluctuation modes acquire a 'mass' k 2 
and decouple. On the other hand, R k J must vanish for large momenta \q\ > k — so 
that the rapid modes remain unaltered and contribute to the functional averages with 
weight exp(— 5fe). Besides, the additional constraints 

lim RV = 0, lim R l i = oo at fixed q (8) 

must be satisfied in order to enforce the correct asymptotic behaviors at the scales 
k = A and k = 0, respectively T k=A ~ S and r fc=0 = T |2"T1 151 12~2"] . 
With the additional term J7J the 'partition functions' 

Zk{],]} = J Vcj)Vi4> exp(-S - AS k + I ' 3<t> + J U) (9) 

become fc-dependent. Finally, the effective T k which is the central object of the NPRG 
procedure is defined as the (modified) Legendre transform of log Z k [j, j] : 

r*fok$ +\ogZ k {j,j] = f - AS,#,# (10) 



Tfe is a functional of the conjugate fields ifi = S\ogZk/Sj and tp — SlogZk/Sj. The 
additional term A<5>fc in Eq. IjlOjl is necessary to set the proper microscopic behavior 
at k = A: T k =A ~ S The RG flow of T k under an infinitesimal change of the 

scale k — or rather s = log(fc/A) — is governed by an exact functional differential 
equation |21II22| (which is derived in the Appendix): 

d s T k = ^Tr f d s R k (f f ] + R k ) ^ . (11) 

In this equation, T k [ip, ip] is the 2x2 matrix of second derivatives of with respect 

to (wrt) ip and ip and [f j?' + Rk}' 1 hence embodies the full (functional) propagator 
associated with the effective theory S k + AS k . 
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Obviously, Eq. cannot be solved exactly and one has to resort to some 

approximations 21 . However, as the approximations used do not rely on the smallness 
of a parameter (see next section), the approach remains non-perturbative in essence. 
In particular, it is not confined to weak-coupling regimes or to the vicinity of critical 
dimensions and is therefore suitable to overcome the limitations of perturbative RG 
schemes. 

4. NPRG for Model A 

To exploit the exact flow equation (|llf) . one has to device an approximation scheme. 
This scheme is based on the construction of an Ansatz for Tk which does not spoil the 
non-perturbative features of the exact equation — and which can be systematically 
enlarged. The formulation of this Ansatz relies on the physics one wishes to probe, 
that is some basic physical insights are necessary. The most common truncation 
consists in expanding Tk in powers of gradients 20 and time derivatives. The accuracy 
and convergence of this approximation scheme have been thoroughly studied in the 
equilibrium context and have shown that quantitatively reliable results can already be 
obtained at the leading order (V 2 ) [21] • For instance, for the three-dimensional Ising 
model, NPRG calculations yield for the critical exponents v = 0.628 and 77 = 0.0443 
at order V 2 E3j , and v = 0.632 and n = 0.033 at order V 4 ]23 : which are in 
close agreement with the 6-loop results v — 0.6304(13) and rj = 0.0335(25) |24|. 
Another useful approximation scheme is the field expansion of Tk- This truncation 
has the advantage of preserving the momentum structure of higher order vertices but 
it approximates the functional structure of the effective potential |2J- The derivative 
expansion is best appropriate for the study of critical physics which is conveyed by 
the large-distance (q — > 0) and long-time (u> — > 0) modes. We hence adopt here this 
approximation scheme and expand Tk at leading order in derivatives — i. e. only terms 
of order V 2 and dt are retained. 

Construction of an Ansatz for Tk at leading order 

The form of the Ansatz for Tk is dictated by the symmetries. Since we consider the 
long-time regime where TRS holds, we want the Ansatz to be invariant under the 
field transformation © (where D is set to unity), which in turn imposes the following 
structure: 

T k ^J)= J otxdt [^X^d^ + ^X^^-X^j^)^ 2 }. (12) 

No higher powers of the response field are allowed at order <9 t due to TRS. Indeed, the 
transformation © would connect a generic 'noise' term tp n X^ n > 2 to higher order 
kinetic terms , : '((9 t -0)™ - -'AT^, j = 0...n which are discarded at order d t . Further 
constraints on the A^'s, i = a,b,c can be deduced from TRS in the same way as 
in Section [3 First one must have A£ = A£ = Xk for the additional contributions 
generated by the transformation © of the 'noise' and the 'kinetic' terms to cancel 
out. As for the linear term in ip, X\. should write as a (field) derivative of a functional 
X^(ip, Vip) = SJ-k/Sip f° r its transform under @ to vanish upon time integration. We 
naturally adopt for Tk the usual equilibrium Ansatz at order V 2 for the Ising model 
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which has been widely studied in the past |^ [231 

F k [4>] = f d<lx {\ Z kO>) [Wf + W)} • (13) 

In this Ansatz, the functional U k embodies the effective potential and the 
renormalization function Z k encompasses the anomalous dimension of the field (see 
below). 

Furthermore, at the scale k = A, T A must identify with the microscopic (bare) 
action © — up to the response field rescaling if) — > ip/Xx — *-e. one has: 

X A = Z A = 1, U A = U{4>) {U{4>) is defined in Eq. ©). (14) 



The Ansatz for Model A at leading order finally writes 

Jx,t 

This Ansatz constitutes the basis of our work 



z k y>) v 2 ^ + i^z fc (^)[v^] 2 



4d^u k (ip)-x k 4> 2 . 

(15) 



Definition of the critical exponents 

We now discuss how the critical exponents can be computed within the NPRG 
approach. In the critical regime, Z k is expected to endow a scaling form Z k ~ k~ nz and 
followingly ip to behave as ip ~ /j(rf-2+»j z )/2 rj^g cr j^ ca i exponent 77 hence corresponds 
to ?7z = —d s \nZk at the critical point [2]. Similarly, Xk is expected to scale as 
Xk ~ fc _J)x at the critical point such that cj ~ k 2 ^ r,z+r,x according to a scaling analysis 
of Eq. (|15l) . Hence, the dynamical exponent z, which by definition characterizes the 
divergence of the time scale following u) — k z 1 is given by z = 2 — r]z + ?7x where 
?7x = -9 s lnX fc . 

Notice that Zk(ip) is a functional of ip whereas the scaling form Zk should be 
a mere (fc-dependent) number. In general, one defines Z k as the value of Zk(ip) at 
a given point ipo, Z k = Zk(4>o). Similarly, though X k here is not a functional, its 
flow equation depends on if), U k {ip), Z k (ip) and their derivatives. Henceforth, the 
notation X k will mean that the corresponding expressions are evaluated for if) = ipo- 
Of course, within the exact renormalization flow, the critical exponents should not in 
fine depend on the choice of if)Q. However any approximation introduces a residual 
dependence and the choice of ipo may become important. The advocated choice (from 
equilibrium studies) is the (running) minimum of the effective potential U k which is 
implicitcly defined by d^U k {ipo) = 0, for it possesses the best 'stability' properties 



Cutoff matrix 

Our last discussion to complete the settings of the NPRG formalism for Model A 
concerns the choice of the cutoff matrix R k . The previous symmetry requirements 
obviously also apply for the quadratic term AS k , which must in particular be invariant 
under TRS (in the stationary regime) . 

The minimal non-perturbative renormalization scheme consists in performing a 
space coarse-graining on the propagator mode ipip, which amounts to considering an 
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off-diagonal cutoff matrix R k with elements 

R 12 = Rf = Q k {q 2 ) and R]}=R 22 = 0. (16) 

This form for Rk is the most natural extension of the equilibrium case — where 
one introduces the scale-dependent quadratic term J ATtk = 1/2 <j>. q Qk(q 2 ) <fiq in the 
(equilibrium) partition function Zk to achieve the splitting of the fluctuation modes 
|21| . We recall that the cutoff function Qk(q 2 ) must decay fastly for large momentum 
modes and behave as fc 2 for slow modes as emphazised in Section [3J A typical cutoff 
function which has been widely used since it allows for analytical results is the 9 cutoff 
introduced by Litim |25| 

Q k (q 2 )^(k 2 -q 2 )9(k 2 -q 2 ). (17) 

It turns out that, even when considering the dynamics, space coarse-graining 
is enough to achieve a proper non-perturbative renormalization program since the 
frequency integrals appear to be convergent and need not be regularized JS]. The 
cutoff matrix (|16fl has hence been adopted in all previous non-equilibrium studies and 
more specifically for reaction-diffusion processes jo]. 

Notice that, on the other hand, one could expect that a time coarse-graining on the 
tpij) propagator may improve the procedure, though it has never been tested. A time 
coarse-graining could be achieved by adding a frequency cutoff R] 2 = £!fc(ia;) on ipip 
modes. But in this case one would have to coarse-grain the noise part correspondingly 
in order to sustain TRS (i.e. the invariance under the transformation © of AiSfc). 
This would amount to introducing an additional cutoff on tpip modes, of the form 
R 22 = — 2i / oj£lk(ioj) ■ The properties of this mixed regularization scheme have never 
been investigated as yet and is left for future work since it represents a great deal of 
numerical efforts. 



5. Flow equations 

The NPRG flow equations for the renormalization functions Uk, Zk and Xk are drawn 
from the exact flow of Tk given by Eq. (|ll|l . According to the Ansatz Ijl5|l . d^Uk can 
be defined by 

o^Uk = -. — tttt hm — 18 

v (2n) d + 1 p,^o §^{p, v) 4>=o V ' 

where the limit of vanishing external momentum and frequency p, v — > means that 
the fields are evaluated in uniform and stationary configurations and the prefactor 
just corresponds to the volume of the system in Fourier space. Similarly, the 
renormalization functions Zk and Xk can be defined by: 

Zk lim d p 2—„ ^ (19) 

( 5 d + 1 (0) p,*-o p 8i>(p,v)5iP(-p,-v) 

X k = gg^ lim div ■ S2T > . (20) 

^+i(0) P ,u^o 5^{ P , v )5^{-p,- v ) ' 

Obviously, X k can alternatively be defined from the noise part as Xk oc 5 2 Tk/Sijj 2 
for a uniform and stationary configuration. One can check that both definitions lead 
to the same flow equation d s Xk which in turn reflects that TRS is preserved by the 
NPRG flow at any scale s. 
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The flow equations of the renormalization functions d^U k , Z k and X k are obtained 
by taking the scale derivative d s of the expressions (|18|l , (|19|l and l|2U[) respectively. It 
is convenient to first rewrite the flow equation d s T k (given by Hill) ) as 

d s T k = i&Trln (ff + , (21) 

where d s (.) = R k S(.)/SR k only acts on the s-dependence of the cutoff elements R% . It 
follows that the field derivatives of d s T k admit simple diagrammatic representations: 

q, (0 

sd s r k i - 



S 2 d s T k 1 ., 



= -d„ 





5ip(p,u)6il)(-p,-i/) 2 

In these graphs, the index 1 or 2 on external legs refers to the corresponding field 
(ip or $ respectively), the n-point vertices correspond to 2 x 2 matrices of rj^ with 
n — 2 (external) fixed indices and two summed over. The propagator lines stand for 
[ri 2) + R k ]-\ which can be easily computed from l|15|) and l|ltj|) . We don't detail 
the full computation of these graphs nor the subsequent explicit frequency integration 
which are lengthy but straigthforward. 

Before giving the resulting flow equations, let us put them in a suitable form 
for the search of fixed point - since we are interested in fine in the scale invariant 
regime. First it is convenient to explicitely express the Ising symmetry by defining 
the functions U k (p) = U k {ip), Z k (p) = Z k (%p) where p = i> 2 j2 is the Z 2 invariant. The 
derivatives of the functions of U k (p) and U k (^p) are simply related: d^U k = ipd p U k , 
d 2 U k = d p U k + 2pd 2 U k . . . and similarly for Z k and Z k . As for the flow equation, 
d s (d p U k ) = l/^d s (d^U k ) for nonzero tfj and d s Z k = d s Z k . Then to absorb any 
explicit dependence on the running scale fc, we introduce the dimensionless quantities 
(according to 1151) ) 

p = k 2 ~ d Z k p 
u(p) = k~ d U k (p) (22) 
z{p) = Z- l Z h {p) 

where the subscript k has been dropped on dimensionless functions. We further 
introduce a dimensionless cutoff function r(y) = Q k (q 2 )/ (Z k q 2 ) where y — q 2 /k 2 
and hence d s Q k = Z k k 2 s(y) with s(y) = —rjzyr{y) — 2 y 2 d y r(y). Finally, the flow 
equations for the dimensionless functions v! and z are given by: 

d s u' = u> {-2 + r) Z ) + (-2 + d + r]z)pu" + 1 (3u" + 2pu"') L d + ^ L 2 + d 

1 

d s z = zr,z + (-2 + d + r i z)pz' + -(z' + 2pz")L d -2pz'{3u" + 2pu"')L d 

\ 2 A/rd 



Ml 



d B In X k 



+ ^{-(l + 2d) pz' 2 L 2 2 +d + 2p{3u" + 2pu'"Y 
+ 4 p z' (3 u" + 2 p u"') M 2+d + 2pz' 2 M 4 4+rf } 
= ^p(3u" + 2pu"'f L d + pz' (3u" + 2pu"') L 2 + d + X -p z' 2 L 4 + d 

(23) 
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where primes denote derivatives wrt p and the so-called thresholds funtions L and M 
are defined by: 



with u" 1 = 2 d n d / 2 T(d/2) and h(y) = y(z+r(y))+u' +2 pu" . By definition (see Section 
SJ, the anomalous dimension rjx = —d s \nXk is obtained by evaluating Eq. (|23|1 at 
the running minimum ipo of the potential Uk, or equivalently at the minimum po of u. 
Similarly, r\z is obtained by solving at the minimum the equation d s z\p a + z' d s po = 
where d s z is given by J22Jl. The additional contribution z' d s pa is generated by the 
running of the minimum implicitely defined as u'(po) = 0. Indeed, the running of u 
implies that its minimum flows according to d s u'(po) — = d s u'\p Q +d s po u"\p which, 
using (|23|) evaluated at p~o yields the expression for d s po. 

We emphasize that, as is to be expected from FDT, the dynamics decouples 
from the statics, i.e. the non-perturbative flow equations (|23[1 for u' and z do not 
depend on Xk (and they identify with those derived in equilibrium with the Ansatz 
11311 21 ). Note that the threshold functions L and M intervening in the equation 
for the anomalous dimensions can be computed analytically upon the choice of the 
cutoff r(y) — (1/y — 1)0(1 — y) (corresponding to i|17|) h which greatly simplifies the 
numerical resolution. The numerical procedure to integrate the flow equations (|23JI is 
detailed in the next section. 

6. Numerical integration of the flow equations 

We will consider different levels of approximation. In a first step — which will be 
referred to as local potential approximation (LPA) — the field dependence of the 
kinetic renormalization function z can be neglected, i. e. only a running coefficient 
is considered. When restoring the p-dependence of z, the corresponding approximation 
will be denoted UZA. 

For both approximations, the numerical procedure to determine the fixed point 
solution of Eq. (|23[) and to compute the critical exponents is the following. We 
sample the field p on a mesh of spacing Ap and discretize the flow equations (|23|l 
using finite differences at order Ap 4 to calculate the p-derivatives of u (and z). 
For the integrals, we either use their analytical expression (whenever available) or 
calculate them numerically using Simpson's rule. We implement an explicit forward 
integration scheme to propagate the solution between scale s and s + As, which turns 
out to be stable for sufficiently small Ap and As. The convergence of the numerical 
procedure when varying Ap and As has been carefully checked. We start out at the 
microscopic scale s = (fc = A) from a quartic bare potential u(p) = X/2(p — pa) 2 
where p\ represents the temperature. We carry through the numerical integration by 
lowering s towards s — ■> — oo. For large bare pa, the system flows to the symmetric 
(high temperature) phase where the (dimensionfull renormalized) minimum of the 
potential pa = k d ~ 2 Z% p~o vanishes, whereas for small p\, it flows to the broken (low 
temperature) phase where po acquires a finite value as s — > — oo. For a fine-tuned initial 
p%, the system is in the critical regime, which corresponds to the effective potential 
u (and z) flowing to a fixed point (scale invariant) form u* (and z*). The critical 




(24) 
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exponents r\ and z can then be computed from the fixed point values of r\z and rjx- 
The critical index v is obtained by linearizing the flow in the vicinity of (u*, z*) and 
determining the (negative) eigenvalue characterizing the unstable (relevant) direction. 
This procedure is carried out within the LPA and UZA approximations. The results 
are gathered in tables and [5] and are commented in the next section. 

7. Results 

The critical exponents for Model A obtained in this work from the NPRG equations 
(|23|l are summarized in Tables ^ an d [21 f° r dimensions d = 3 and d = 2 respectively, 
and compared with results ensuing from other field theoretical methods and Monte 
Carlo simulations. 

Let us first comment on the equilibrium exponents r\ and v. As emphasized 
in Section [5] the dynamics decouples from the statics in Eqs. (|23[) — as expected 
from TRS in the stationary regime. As a consequence, the equilibrium exponents 
computed here should match (up to numerical accuracy) those obtained in earlier 
NPRG works on the (equilibrium) Ising model. The three-dimensional Ising model 
has been thoroughly investigated within the NPRG framework as a testing ground of 
the method |2f>| . In particular, critical exponents have been computed using the cutoff 
function (fTTfl in j^Hl 123 (though with different numerical procedures) . The exponents 
v and 77 we obtain in d = 3 precisely reproduce these values both at LPA and UZA. 

We know from these previous studies that the accuracy can be improved by 
optimizing the choice of the cutoff function. At order V 2 in derivatives, optimized 
exponents are v = 0.6281 and rj = 0.0443 [2E]j which are already in close agreement 
with the 6-loop calculations |2U or Monte Carlo simulations [35]. However, since the 
determination of r\ is related to the momemtum structure of the two-point correlation 
function, its accuracy is poorer than that of v. A better accuracy on r\ requires 
to compute the next order V 4 in derivatives which yields 77 = 0.033 23 . In two 
dimensions, far fewer NPRG results are available. A calculation with cutoff function 
Ijl7|) has been achieved in and both results are in close agreement. As for in d = 3, 
the determination of i] remains poorer than that for v at order V 2 in derivatives. 
However, the two-dimensional Ising model has not been systematically investigated 
and neither optimized nor order V 4 exponents have been determined in d = 2. 

We can now come to the new part of this work which concerns the dynamics. 
The situation for z is very different from that of the equilibrium critical exponents. 
For the dynamics, no high-loop expansions or exact results in d = 2 are available. 
Furthermore, results from MC simulations appear to be rather scattered especially in 
d = 2. The values reported in Tables and [2] seem to be accepted as reference values 
|T%) . On the field-theoretical side, the determination of z is very sensitive to the choice 
of the resummation scheme since only a few orders are known. Various resummation 
schemes have been studied and we quote here the latest results obtained in |50|. 

Our results are in reasonable agreement with these various estimates. This is 
one of the key point of the NPRG approach and a central motivation for this work: 
the leading order in derivatives appears to already provide a reliable determination of 
physical quantities, as outlined above. Notice that the variation on z between LPA and 
UZA is not meaningful and merely provides an estimate of the error. The reason is that 
going from LPA to UZA does not amount to enriching the Ansatz for the dynamical 
part since in both cases only a running coefficient Xk is allowed by TRS. Hence even at 
UZA, the accuracy on z remains poorer than that for the equilibrium exponents and 
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d = 3 


V 


V 


z 


Ref. 


LPA 


0.65 


0.11 


2.05 


[this work] 


UZA 


0.63 


0.05 


2.14 


[this work] 


FT 


0.6304(13) 


0.0335(25) 




Eg 


MC 


0.6297(5) 


0.0362(8) 




Em 


FT 






2.0237(55) 


Em 


MC 






2.032(4) 2.055(10) 


mums 



Table 1. Critical exponents of Model A in d = 3 from the different NPRG 
approximations (LPA and UZA) computed in this work, compared with results 
from field theoretical methods (FT) and Monte Carlo simulations (MC). 



d = 2 


V 


V 


z 




LPA 


0.78 


0.43 


2.15 


[this work] 


UZA 


1.1 


0.37 


2.17 


[this work] 


exact 


1 


1/4 






FT 






2.0842(39) 


Em 


MC 






2.1667(5) 





Table 2. Critical exponents of Model A in d = 2 from the different NPRG 
approximations (LPA and UZA) computed in this work, compared with results 
from field theoretical methods (FT) and Monte Carlo simulations (MC). 



the rapidity of convergence on z can not be tested within these approximations. One 
would need to implement the next order in time derivatives to improve the accuracy 
on z, which is rather costly. Alternatively, one could modify the regularization scheme 
and resort to a frequency and noise coarse-graining as mentioned in Section which 
is likely to yield better results, but is yet to be investigated. 



8. Summary 

In this work, we have studied the critical dynamics of Model A within the NPRG 
formalism. We have in particular detailed the device of an appropriate approximation 
scheme preserving the symmetries, the derivation of the NPRG flow equations and 
their resolution. Using a very simple Ansatz, that is at the leading order in derivatives, 
we have obtained a reliable estimate for the dynamical exponent z: z — 2.09(4) in d = 3 
and z — 2.16(1) in d — 2. The fact that the leading order already yields quantitative 
results is a generic feature of the NPRG approach which makes it particularly 
powerful. This feature is the central motivation for the emphasis put throughout 
this work on the methodological part: a leading order NPRG calculation can already 
enable one to investigate nontrivial problems. Restricting to non-equilibrium critical 
phenomena, NPRG studies have indeed brought out new non-perturbative properties 
of reaction-diffusion processes |10| and allowed to tackle interface growth problems 
|34| . We hence believe the approach to be useful to investigate many other non- 
equilibrium scaling phenomena. Of course, as the application of NPRG techniques 
to non-equilibrium statistical physics and dynamics is very recent, a great deal of 
systematic studies remain to be done as well to test the efficiency of the different 
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Ansatz and regularisation schemes out-of-equilibrium, which will be the goal of future 
works. 

The authors are indebted to B. Delamotte for fruitful exchanges all through this 
work. The authors also wish to thank G. Biroli, A. Lefevre and N. Wschebor for 
enlightening discussions. 



9. Appendix: NPRG flow equation for T k 



This appendix is dedicated to the derivation of the exact flow equation Ullfl for 
T k . (A similar derivation can be found for equilibrium systems in 21 ). In this 
appendix, we use the shorthand x = (x,t), vectors are denoted by capital letters 
(eg. = [ij),i{>], J = \j, j]) and for functional derivatives we introduce the notation 
Sj x = [5/5j(x), S/5j(x)] (and similarly for <5* x ). Overhead hat symbols are used for 
(2 x 2) matrices. 

The variation of T k (at fixed 'J) under an infinitesimal change of the running 
scale k follows from Eq. i|l(J|) : 



*(x).9 fc i? fc (x,x'). t *(x') 



dkTk^k]^ = -d k W k [J]\ 9 + / 9 fct 7(x).^ fe (x) 



(25) 

where W k = \nZ k . The variation with k of W k at fixed ^ is related to that of W k at 
fixed J by: 

■% y W k . (26) 
The expression of <9fc>Vfc| „ is obtained by taking the derivative of © wrt k: 



d k W k 


= d k W k 






* 





d k In Zi 















J 







-S k - AS k +' <$>3 



{-\J 6 Jx .d k R k (x,*') 
n. iv, 



= d k W k e 



After expressing the derivatives 5j and t 5j of exp (W k ) and dividing out by exp (W k ) 
one obtains: 



d k W k 



J 



I J {^W fc .5 feJ R fe (x,x'). t ^W fe +Tr [d k R k (x,x.') .5 Jx (% x ,W k )]} . 



. (27) 

The last term in the right hand side is the matrix of second (functional) derivatives of 
Wfe which we denote (x, x'). The flow equation of T k follows from inserting (|26[l 
and I23 in which yields 



d k T k m\^ = ^Tr I 4iJ,(x,x').wf ) (x,x'). 

* Jx.x' 

This equation can be conveniently expressed in a closed form upon inverting 



(28) 



(2) 



; . The inverse of Wl z> can be obtained by taking a functional derivative t S^ x , of 



the definition *(x) = 5 Jx W k : 

T.^ +1 )(x,x')= / % x ,J(y).wi 2) (y,x). 

J y 



(29) 
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The matrix '<!>* , J{y) is simply given by two successive derivatives 6y and *<!>* , of 
Eq i fTUj l. which yields : 



t ^J(y)=ri 2) (x',y)+i? fc (x',y), (30) 



where Y k (x, x') denotes the matrix of second (functional) derivatives of IV Finally, 
inserting 13U|) in the flow equation (|28|l yields the advocated equation: 



which can be Fourier transformed and underlies the NPRG calculations of this work. 
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